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ABSTRACT 

The average optical-UV continuum shape of QSOs has been investigated using spec- 
tra of 62 QSOs having good relative photometric calibrations. The QSO spectra were ex- 
tracted from two complete color selected samples in the magnitude intervals B m 16 — 20. 
The analysis was performed fitting power-law continua f v oc v a in well defined rest-frame 
wavelength intervals after removing regions of the spectrum affected by strong emission 
lines or weak emission bumps. The average slope in the rest-frame optical-UV region 
1200 — 5500 A shows a rapid change around the 3000 A emission bump with a ~ 0.15 
longward of it and a ~ —0.65 at shorter wavelengths. Although these average slopes 
have been obtained using spectra of QSOs with different luminosities and redshifts, there 
are no significant correlations of the average spectral index with these quantities. For 
a few QSOs in the sample we were able to measure the same softening of the spectral 
shape within the individual spectrum. These results have significant consequences on 
the estimate of the cosmological evolution of the optically selected QSOs as they affect, 
for instance, the fc-corrections. New /c-corrections in the B, V, R and Gr bands were 
computed. The derived average spectral shape in the optical-UV band puts interesting 
constraints on the expected emission mechanisms. 
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1. Introduction 

The study of the shape of the QSO continuum has 
a great importance both for the identification of the 
emission mechanisms responsible for the observed ra- 
diation and for the estimate of the QSO cosmological 
evolution. 

' Lhe continuum spectrum is conventionally param- 



ete 'izcd in the optical- U V region by a power law of 



the type /„ oc is". This is a local approximation: in 
the recent years the availability of IR, UV and X-ray 
data has allowed to measure the energy distribution 
of QSOs in a wide spectral range, showing large devi- 
ations from a single power-law approximation. Even 
in the optical-UV region, composite spectra derived 



from complete color selected samples (Cristiani & Vio 



1990, Francis et al. 1991) show significant departures 
from a single power-law. 

Another remarkable occurrence is the variety of 
results obtained in the many studies of the shape 
of the optical-UV continuum spectrum, in particu- 
lar about the estimated average optical-UV spectral 
index, a, and its dispersion. Roughly the intrin- 
sic optical-UV spectral index ranges from -1 to 0. 
O'Brien et al. (1988) and Sargent et al. (1989) ob- 
tain —0.7 < a < —0.6 on two large samples; Oke 
& Korycansky (1982), Oke et al. (1984), Baldwin 
et al. (1989) and Cheng et al. (1991), found values 
in the range —1 < a < —0.5 from smaller samples. 



dex 



~ —0.3 (Banders et 


al. 198£, Francis et al. 


Webster et al. 1995 




Francis 1996|), in agree- 



ment with the expected slope provided by free-free 



emission models (Barvainis 1993) 



However, most of the differences among the results 
quoted in the literature might be due to the different 
spectral intervals used to estimate the average con- 
tinuum slope. 

No substantial agreement exists also about the dis- 
persion of the observed spectral index, er Q , and its 
origin: obscuring dust (Webster et al. 1995), intrin- 



sic differences in the spectral energy distribution from 
QSO to Q SO ( [Sargent et al. 1989| , [Serjeant fc Rawl- 



ngs 1996), or simply measurement errors (Cheng et 
il. 199l|). Some authors claim no intrinsic disper- 



sion (Banders et al. 1989] ), others state a large one 
(Elvis et al. 1994): so far the observed values of the 
dispersion range from 0.2 to 0.5. 

Determining the average shape of the optical-UV 
emission and whether it is intrinsically the same for 



all the QSOs is not a worthless question. The choice 
of a = —0.7 rather than -0.3 to calculate the k- 
corrections leads to increase the luminosities of QSOs 
at z ~ 2 of a factor w 2. An appreciable disper- 
sion (<t = 0.3), as well, introduces uncertainties in 
the luminosity and biases the estimate of the cosmo- 
logical evolution of the luminosi ty function (piallongo 



1997) 



& Vagnetti 1992, Francis 1993, La Franca & Cristiani 



In this paper we present a spectroscopic study of 
a set of 62 objects extracted from a complete sample 
of radio-quiet QSOs. We determine the values of the 
mean spectral index in four different ranges of rest 
frame wavelength, with a particular attention to the 
correlations between a and the redshift or the lumi- 
nosity. 

We describe in Sect. 2 the data set and in Sect. 3 
the method we used to perform the power law fit to 
the continuum. Finally, in Sect. 4, we show the ob- 
served dependence of the spectral index on the mean 
rest frame wavelength and present new fc-corrections 
to estimate the evolution of the QSO luminosity func- 
tion (LF). 

2. Data sample 



The QSOs studied in the present work belong to 
the Selected Area 94 (La Franca et al. 1992; < 
19.9) and the Homogeneous Bright Quasar Survey 
(HBQS, |Cristiani et al. 1995| ; 15 < m B <18.75). 
Both samples identified quasar candidates on the ba- 
sis of UV criteria (complementary objective-prism 
data have been also used in a few subareas to cross 
check the effectiveness of the selection). From this 
sample a subset of 62 QSOs was extracted with a 
good relative photometric calibration and with a lu- 
minosity and redshift distribution similar to the orig- 
inal set (Figure 1). It is to notice that the spectral 
properties of our sample are representative of QSOs 
selected with an UV excess and a bias against radio- 
loud quasars showing very red colors can be present 
( |Webster et al. 1995| , [Francis 1996| ). On the other 
hand it is not known at present if radio-quiet dust- 
reddened objects represent a considerable fraction of 
the QSO population. 

The observations of the QSOs have been carried 
out at the 3.6m, 2.2m ESO/MPI and 1.5m ESO tele- 
scopes at La Silla, equipped either with the ESO 
Faint Object Spectrograph and Camera (EFOSC1 - 
EFOSC2), or with the Boiler and Chivens Spectro- 
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Fig. 1. — Luminosity redshift plot of the QSO 
sub-sample studied. The apparent magnitude loci at 
tub = 15 and tub = 20 are shown (Jiq = 50; qo = 0.5). 



graphs. The slit had typically a width at least 1.5 
times the seeing FWHM and in general was placed 
at the parallactic angle. The detectors were always 
CCDs. The resolution of the spectra ranged between 
10 and 30 A. The reduction process used the stan- 



dard MIDAS facilities (Banse et al. 1983) available 
at the Padova Department of Astronomy and at ESO 
Garching. 

The raw data were sky-subtracted and corrected 
for pixel-to-pixel sensitivity variations by division by 
a suitably normalized exposure of the spectrum of an 
incandescent source. Wavelength calibration was car- 
ried out by comparison with exposures of He-Ar, He, 
Ar and Ne lamps. Relative flux calibration was finally 
achieved by observations of standard stars listed by 
Dke 1974 and Stone 1977. The airmass of the obser- 



vation was always below 1.6 and has been corrected 
on the basis of the standard La Silla extinction tables. 

3. Fitting the continuum 

As discussed in the previous sections the shape 
of the continuum component of the QSO spectra 
changes considerably over different frequency bands. 

In the optical-UV band the continuum is usually 
parametrized as a power law: /„ cx v a . Some results 
dCristiani fc Vio 1990| , [Francis et al. 199l| ) show that, 
even in the optical-UV range (1200-10000 A), the con- 
tinuum spectrum of individual QSOs is not well rep- 
resented by a single power law (i.e. a = constant over 
the whole band). However these results are obtained 



with the technique of the composite spectra, based on 
the simplifying assumption that each spectrum is rep- 
resentative of the whole QSO "class" and an ensemble 
mean can be meaningfully applied. Such an hypothe- 
sis may result too simplistic and even misleading: it is 
well known, for example, that the average of a num- 
ber of power-laws of different spectral index is not a 
single power-law. 

In this study we have estimated the local power-law 
shape of the individual QSOs in four ranges, namely 
the rest-frame intervals: 1400 - 2200 A, 2150 - 3200 
A, 2950 - 4300 A, 3900 - 5500 A. Some spectra of the 
extracted set of QSOs cover a large frequency range, 
so in few cases we have, for the same QSO, two dif- 
ferent rest-frame wavelength intervals. We have ob- 
tained 78 spectral ranges from the extracted set of 62 
QSOs. 

For each spectral range the regions certainly con- 
taminated by strong emission were eliminated and 
an automatic algorithm was applied to estimate the 
power law index a over the spectral regions described 
above. The number of available regions in individ- 
ual QSO spectra depends on the redshift and on the 
instrumental response of the spectrographs used. 

The method consists in a first selection of the prob- 
able continuum window regions within each spectral 
interval. These continuum windows are selected on 
the basis of the composite emission line spectra pub- 
lished by Cristiani & Vio (1990) and Francis et al. 
(1991). 

The local power-law continuum is determined by 
a linear regression on the local minima of the spec- 
trum selected after appropriate smoothing to improve 
the S/N ratio. Each minimum is weighted by the 
reciprocal of the local second derivative of the spec- 
trum, to minimize the effects of occasional absorp- 
tions. The weight is doubled for those local minima 
which, following the results on the composite spectra, 
are within identified continuum windows (see Table 
!)• 

4. Results 

4.1. Continuum slope 

The mean value of the spectral index, SIT, aver- 
aged over the whole sample, is -0.33 with a root mean 
square dispersion of 0.59, in agreement with the value 
found by Francis (1993) for the Large Bright Quasar 
Survey (LBQS). This value of the dispersion, however, 
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Fig. 2. — Examples of QSOs' spectra with fitted continuum and resulting spectral index in the given rest wavelength 
interval. 
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although supported by Elvis et al. (1994), is signif- 
icantly higher than that (« 0.3) found by Sanders 
et al. (1989) and recently by Francis (1996). If the 
power-law index is a function of the wavelength, the 
different results can be understood in terms of mea- 
surements in different rest-frame spectral regions: in- 
deed, under the hypothesis of significant deviations 
from a single power law, we do expect a large disper- 
sion over the whole sample. 

The mean values of the continuum slopes obtained 
for each of the four rest wavelength intervals described 
above are reported in Table 1 together with their dis- 
persions. For 7 QSOs it was possible to fit the contin- 
uum only in the range 1700 — 2700 A. Figure 3 shows 
the dependence of the spectral index from the mean 
rest wavelength of the fitted spectrum: noticeable the 
presence of a break in the values from cE^T ~ —0.65 to 
5^7 ~ 0.15 at A ~ 3200 A. In this case the disper- 
sions are lower (see Table 1) than found for the whole 
sample. 



-2 
1000 



3000 4000 
neon rest wavelength (Angstrorr 



Fig. 3. — Dependence of the mean spectral index 
from the mean rest wavelength of the fitted spec- 
trum. Along the y-axis are plotted the errors in a, 
while along the x-axis the average baseline of the fit 
is shown for each rest wavelength interval. 

To represent the overall shape of the average QSO 
continuum as a gradual change of the local slope as a 
function of wavelength we have assumed a linear trend 
of a v vs X rest (Figure 3), i.e. = a(j^) + b (where 
the wavelength is expressed in Angstrom) ; we derived 
a = 1.5 ± 0.3 and b = —1.5 ± 0.2 by a linear fit with 



errors in both coordinates (see Fasano & Vio 1988). 
In this case, because of a v — d ^ g J , the optical- 
UV continuum is expressed by a more complicated 



function of the wavelength: 



LogF u 



vLogv 



Logv 



(1) 



where k 



(c is the speed of light). 



4400 inlO 

No significant correlation between a v and the red- 
shift or the luminosity was found (Table 2), in agree- 
ment with the results obtained by Richstone & Schmidt 
(1980), Oke et al. (1984) and Cheng et al. (1991). 
In all the four wavelength ranges we have 



0.4 ± 0.4 and 



da. u 
dM B 



dz 



< 



|< 0.1 ± 0.1, with a probabil- 



ity greater than 15% of no correlation between a v 
and redshift or luminosity. Since has been found 
roughly constant below and above 3200 A, we stud- 
ied the correlation with z and Mb in the two larger 
wavelength intervals 1400 - 3200 A and 3200 - 5500 
A (Figure 4) , combining the values of a found in each 
of the intervals reported in Table 1: even in this case 
no significant correlation was found. 

To avoid spurious results due to the luminosity- 
redshift correlation present, to some degree, in flux- 
limited QSO samples, we repeated the correlation 
analysis in restricted intervals of luminosity and red- 
shift. The mean values of a v and their degree of cor- 
relation with z and Mb do not change substantially 
if we choose objects of the sample having luminosities 
or redshifts in the following intervals: —27.5 < Mb < 
—25.5, 0.5 < z < 1.1. These results suggest that the 
variations of the spectral index with A res t are intrin- 
sic and do not depend on the distribution of QSOs in 
the redshift-luminosity plane. 

As a further check for thirteen objects it was pos- 
sible to measure the change in the average continuum 
slope around 3000 A by fitting the individual continua 
over two different (contiguous) wavelength baselines: 
five QSOs over AA2150 - 3200 A and AA2950 - 4300 
A, eight over AA2950 - 4300 A and AA3900 - 5550 A. 
The resulting mean difference in the spectral indices 
between the two spectral intervals, averaged over each 
of the two sets of QSOs, are, in both cases, close to the 
difference in the values of the mean spectral indices 
found for the whole sample. 
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rest wavelength interval 


continuum windows 


OL v 


a 


number of QSOs 


1400-2200 A 


1420-1480 A 
2150-2300 A 


-0.61 


0.37 


19 


1700-2700 A 


2150-2300 A 
2600-2700 A 


-0.89 


0.32 


7 




2150-2300 A 








2150-3200 A 


2600-2700 A 
2950-3100 A 


-0.70 


0.35 


19 


2950-4300 A 


2950-3100 A 
3900-4300 A 


0.23 


0.48 


22 


3900-5500 A 


3900-4300 A 
5400-5500 A 


0.07 


0.39 


11 



Table 1 : Mean values and dispersions of the continuum slopes obtained for each of the four rest wavelength intervals 
in which the variations of the spectral index a are within the errors of the fit; also the result for the interval 1700 — 
2700 A (7 QSOs) is shown. In the second column are reported the continuum window regions (identified following 
the results on the composite spectra) falling in each rest wavelength interval. 
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rest wavelength interval 


1400-2200 A 


1700-2700 A 


2150-3200 A 


2950-4300 A 3900-5500A 


dz 


0.43 ±0.36 




-0.11 ± 0.27 


0.42 ±0.47 0.30 ±0.63 




0.277 




0.101 


0.199 0.158 


P(r a ,z) 


0.25 




0.68 


0.37 0.64 


da u 
dU B 


0.07 ±0.07 




-0.02 ±0.06 


-0.07 ±0.08 -0.06 ±0.11 


r a .M B 


0.246 




0.069 


0.181 0.171 


P{r a .M B ) 


0.31 




0.78 


0.42 0.61 


dz 




0.11 ±0.12 




0.45 ±0.33 


Ta,z 




0.146 




0.239 


P(r a ,z) 




0.34 




0.18 


da u 
dM B 




-0.00 ±0.04 




-0.07 ±0.06 


r a,M B 




0.013 




0.187 


P{r a ,M B ) 




0.93 




0.30 



Table 2: Correlation between the spectral index a v and both the redshift and the luminosity; r is the correlation 
coefficient and P(r) is the probability of no correlation. The absolute magnitude was calculated assuming ha — 50 
and go = 0.5. 
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Fig. 4. — Correlation between the spectral index and both the redshift (left) and the luminosity (right). Different 
symbols indicate values of a belonging to different rest-frame intervals: 1400 — 2200 (diamond); 1700 — 2700 
(square); 2150 - 3200 (cross); 2950 - 4300 (circles); 3900 - 5500 (triangles). Here is assumed ho = 50 and 
go = 0.5. 



4.2. fc-corrections and the evolution of the 
QSO luminosity function 

The shape of the continuum spectrum has an im- 
portant role in determining the fc-corrections and 
therefore the luminosity function of the QSOs. In 
fact, because of the different redshifts of the observed 
objects, to evaluate properly the QSO monochromatic 
luminosity function and its evolution, the observed 
flux of the QSOs has to be corrected by a quantity, 
the fc-correction, that depends on the shape of the 
spectrum. 

The fc-correction is given by the usual formula 



k = 2.5log(l + z) + 2.5/05 



f °°F(\) S (\)d\ 
J^F(^- z )s(\)dX 



(2) 



if-corrections in the Johnson's B, V, R and Gunn's 
Gr bands (Table 3) were calculated using equation (2) 
with a composite spectrum (Figure 7) built as follows: 
we computed the continuum spectrum from equation 
(1) and the emission component averaging the residu- 
als, respect to their own continuum, of each spectrum 
of the sample. The fc-corrections in the blue bandpass 
were compared (Figure 5) with the ones derived from 
two composite spectra having a power law continuum 
with constant spectral indexes, respectively a = —0.3 
and a = —0.5 and with the fc-corrections given by 
Cristiani & Vio (1990). 
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Table 3: X-corrections in the B, V, R (Johnson) and Gr (Gunn) bandpasses as derived from our composite 
spectrum. 



9 



It is to notice that for z < 1 there is a good agree- 
ment between our fc-corrections and the ones derived 
from a spectrum with a constant slope a v = —0.3, 
whereas the difference rises to about 0.2 magnitudes 
at z ~ 2. This suggests that the hypothesis of a flat 
single-power-law spectrum leads to an underestimate 
of the luminosity of the QSO at high redshifts. Our 
results are in better agreement with the fc-correction 
computed by Cristiani & Vio (1990) although in this 
case lower luminosities are predicted at z < 2 by our 
curve. 




Fig. 5. — if -corrections in the B bandpass calculated 
in the case of our synthetic spectrum (solid) and spec- 
tra having single power-law continuum with spectral 
indexes -0.3 (dotted) and -0.5 (dot-dashed) compared 
with fc-corrections given by Cristiani & Vio (1990) 
( dashed) . 

This new fc-correction was used to estimate, by a 
maximum likelihood method, the optical luminosity 
function $(Mb, z) for a large sample of ~ 1000 QSOs, 
taken from various surveys in the literature as de- 
scribed in La Franca & Cristiani (1997). A double 
power law luminosity function and a pure luminosty 



evolution (PLE) model of the type (see La Franca & 
Cristiani 1997 for details): 
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were assumed, where a and /3 are the slopes of 
the faint-end and bright-end luminosity function, re- 
spectively. $(Mb,z) is in units of Mpc _3 mag _1 . 
In this model the evolution is specified by the red- 
shift dependence of the break magnitude (M B (z) = 
M* B - 2.5klog(l + z)), where M* B = M£(z=0). This 
relationship corresponds to a power law evolution of 
the luminosity: 



L*(z) = L*(l + z) K 



(4) 



The results are shown in Table 4 and Figure 6. The 
luminosity function computed using our fc-correction 
has a value of the faint-end slope a between the values 
of the faint-end slope derived from the single power- 
law spectra, while the bright-end slope (3 is the same 
(within the errors). Moreover, the z = break lumi- 
nosity (M B ~ —22.5) does not change substantially 
with respect to the LF computed with the assump- 
tion of a single power-law slope, but the parameter 
k, which gives the rate of the cosmological evolution, 
goes from the value 2.90 for a spectral index -0.3 to 
3.12 for our spectrum, and the probability of the fit 
improves noticeably (Table 4). 

Comparing our luminosity function with that com- 
puted by La Franca & Cristiani (1997) using the Cris- 
tiani & Vio (1990) fc-correction we confirm the flatten- 
ing at low redshift of the bright part of the luminosity 
funtion in comparison with a Pure Luminosity Evolu- 
tion model. This is indicated by the poor fitting prob- 
abilities given by the two dimensional Kolmogorov- 
Smirnov test in the range (0.3 < z < 0.6) (see Ta- 
ble 4). However, the larger differences between our 
and Cristiani & Vio (1990) fc-correction are at z ~ 2, 
where the discrepancy is of 0.2 magnitudes. This 
is at the origin of a evolutionary parameter which 
is about 2cr lower, resulting in a volume density of 
QSOs at z ~ 2 which is lower by a factor about 1.5 
at M B 27.5. 

Our more accurate estimates of the dispersions of 
the average spectral index as a function of wavelength 
(a = 0.3 — 0.4 from Table 1) can also be included 
in the evaluation of the evolution of the luminosity 
function. Uncertainties in the estimate of the average 
spectral slopes are known to produce a smaller cos- 



mological evolution of the luminosity function ( Gial- 
longo fc Vagnetti 1992] ). In this case a value for the 



evolutionary parameter as small as k = 2.94 has been 
found (see Table 4). The probability of the PLE evo- 
lutionary models with the new fc-corrections are in 
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lcr errors 




±0.07 


±0.13 


±0.07 
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Table 4: Parameters for luminosity function with different continuum shapes: (a) PLE model (qo = 0.5); (b) PLE 
model including a dispersion in the distribution of the spectral slope a = 0.4 (go = 0.5). The last two columns are 
respectively the % 2 test probability and the two-dimensional Kolmogorov-Smirnov test probability (see details in 
La Franca & Cristiani, 1997). 



general higher than previously found (P(x 2 ) ~ 0.5), 
showing how the selection of a more accurate average 
spectral shape can discriminate between different evo- 
lutionary models for the QSO luminosity function. In 
particular, there is no evidence for any redshift cutoff 
in the redshift evolution of the luminosity function as 
introduced by Boyle (1992) at z = 1.9. 
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Fig. 6. — The QSO luminosity function. The points 
connected with a continuous line represent the obser- 
vations (see La Franca & Cristiani, 1997 for details); 
the dotted lines are the best-fitting Pure Luminosity 
Evolution model. Error bars correspond to la confi- 
dence intervals. 



4.3. Modelling the emission mechanism 

As an example of the consequences that more ex- 
tended works of this type can have in constraining 
thermal emission models for the QSO continuum, we 
have compared with thermal models based on the 
presence of an accretion disk around a massive Ker r 
black hole (10 6 - 10 9 M Q ) (|Czerny fc Elvis 1987|) . 



These models have been already used to fit the con- 
tinua of QSOs in a wide wavelength in terval extending 
to the X-ray band ( Fiore et al. 1995 ). 

The emission comes from a geometrically thin ac- 
cretion disk around a massive black hole with accre- 
tion rates m ranging in the sub-Eddington regime. 
For accretion rates much higher than the Eddington 
one the disks become slim and additional cooling pro- 
cesses should be included. In this illustrative example 
only sub-Eddington thin disks are considered. 

The average shape has been normalized to Mb ~ 
—24 since higher luminosities require emission in a 
super-Eddington regime which is not allowed by the 
model. A qualitative comparison seems to indi- 
cate that only models with a large black hole mass 
(M ~ 10 9 ' 5 M Q ), an average inclination along the 
line of sight fi = cosO = 0.5 with an interme- 
diate accretion rates m = 0.3, or smaller inclina- 
tion (/i = cosO = 0.75) and higher accretion rate 
(m = 0.8) can reproduce a change of slope Act w 1 
at Aw 3200 A (Figure 7). Models with smaller black 
hole masses (M < 10 7 ' 5 Mq) provide much lower lu- 
minosities, while models with masses between 1O 8 M0 
and 1O 8 5 M0 present a softer continuum shape than 
the observed one. For black hole masses M m 10 9 — 
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Fig. 7. — Composite spectrum with the continuum 
component (solid), as estimated in the present work, 
compared to the synthetic continuum spectrum de- 
rived from thermal emission models of an accretion 
disk around a massive Kerr's black hole. The best fit- 
ting model (M = 1O 9 5 M , m = 0.8 and cos9 = 0.75; 
see Fiore et al., 1995, for details) is shown (dashed). 

10 9 5 M Q , high accretion rates, and smallest inclina- 
tion angles (cos9 — 1 corresponds to a face-on disk) 
there is a good agreement for Logv < 15.2, but the 
large masses involved in the accretion process pro- 
duce a drop of the luminosity in the blue part of the 
spectrum. 

5. Conclusion 

In this paper we have analysed the spectra of 62 
QSOs in the magnitude interval B = 15 — 20 having 
good relative photometric calibrations. The analysis 
was performed fitting power-law continua /„ oc v a 
in well defined rest-frame wavelength intervals after 
removing regions of the spectrum affected by strong 
emission lines or weak emission bumps. 

The main results of the paper can be summarized 
as follows: 

• The mean values of the continuum slopes ob- 
tained in the spectra of 62 QSOs change from 

a ~ 0.15 at A > 3000 A to a 0.65 at 

A < 3000 A. 

• Assuming a linear change of the local slope with 
the rest-frame wavelength we obtain the follow- 
ing relation a = 1.5(A/4400) - 1.5 



• There is no correlation in the sample between a 
and redshift and/or luminosity. 

• New fc-corrections, which include the contribu- 
tion of emission features, are derived for the B, 
V, R (Johnson) and Gr (Gunn) bands. 

• The new fc-corrections are used to estimate the 
evolution of the QSO luminosity function. As- 
suming PLE models, the double power-law lu- 
minosity function evolves as L oc (1 + z) 31 or as 
low as L oc (1 + z) 2 - 9 if we take into account the 
dispersion in the distribution of spectral shapes. 

• A qualitative comparison of the QSO average 
continuum shape with thermal models based on 
the presence of a thin accretion disk around a 
massive Kerr black hole seems to indicate that 
only models with a large black hole mass (M ~ 
1O 9 5 M ) can reproduce a change of slope Aa w 
1 at A « 3200 A. 
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